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Abstract 

We study the corrections to adiabatic dynamics of two coupled quantum 
dot spin-qubits, each dot singly occupied with an electron, in the context of 
a quantum computing operation. Tunneling causes double occupancy at the 
conclusion of an operation and constitutes a processing error. We model the 
gate operation with an effective two-level system, where non-adiabatic transi- 
tions correspond to double occupancy. The model is integrable and possesses 
three independent parameters. We confirm the accuracy of Dykhne's formula, 
^ ■ a nonperturbative estimate of transitions, and discuss physically intuitive con- 

ditions for its validity. Our semiclassical results are in excellent agreement with 
numerical simulations of the exact time evolution. A similar approach applies 
to two-level systems in different contexts. 
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1 Introduction 

Quantum information processing is an active and fascinating direction of research with 
participation from various fields of physics and neighboring scientific disciplines pQ. 
This extraordinary interest has generated a fairly vast amount of theoretical and ex- 
perimental studies. Possible experimental realizations of quantum information pro- 
cessing are presently being investigated. Among the different approaches, those in a 
solid state setting are attractive, because they offer the potential of scalability — the 
integration of a large number of quantum gates into a quantum computer once the 
individual gates and qubits are established. With that in mind, several proposals for 
using electron and/or nuclear spins in solid state systems have been put forward in 
recent years j3 El HI HI EH U\ ■ Specifically, in Ref . [2] it was proposed to use the spin of 
electrons residing in semiconductor quantum dots as qubits |H} li^ lTm ITU IT^ITBIITH IT^]. 
In this paper we revisit the quantum dynamics of gate operations between qubits of 
this type. Such two-qubit operations are performed by varying the amplitude of elec- 
tron tunneling between the dots via external electric potentials. In a generic scenario, 
the tunneling amplitude between the dots is zero (or, more precisely, exponentially 
small) before and after the gate operation, while it is finite and appreciable during 
such a process. Thus, the typical time dependence of the tunneling amplitude is a 
pulse roughly characterized by its duration, amplitude, and ramp time (see Figure d) 
During such a pulse, the tunneling amplitude is finite and essentially constant, and 
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both electrons can explore the total two quantum dot system. Therefore, their indis- 
tinguishable fermionic character is of relevance fDJEHEl- m particular, in such gate 
operations entanglement-like quantum correlations arise which require a description 
different from the usual entanglement between distinguishable parties (Alice, Bob, ...) 
in bipartite (or multipartite) systems. In such a case the proper statistics of the in- 
distinguishable particles has to be taken into account [TQl EE HZ] • 

Another important aspect of having a finite (as opposed to infinitely high) tunnel- 
ing barrier between the dots is that it necessarily leads to (partially) doubly occupied 
states in the two-electron wave function, i.e. contributions to the wave function where 
both electrons are on the same dot (having different spin) occur with finite ampli- 
tude. Doubly occupied states which arise as the result of a measurement after the 
gate operation, destroy the information in those qubits and lead to errors in the infor- 
mation processing. Therefore, it is desirable to reduce the probability of such errors, 
i.e. the occurrence of doubly occupied states, in the resulting two-electron state af- 
ter the gate operation, while it is necessarily finite during the operation j§] EI] . If 
the error probability can be sufficiently reduced, error events can be tolerable and 
handled with quantum error correction schemes. An effective way of guaranteeing 
error suppression is to maintain nearly adiabatic time evolution. Doubly occupied 
states then correspond to corrections to adiabatic evolution, which are often called 
"non-adiabatic transitions". Numerical simulations [TU] have shown that the adia- 
batic region, in terms of the pulse parameters such as ramp time and amplitude, is 
rather large. On a heuristic level, this numerical result is plausible on the basis of 
the classic papers on adiabatic quantum motion in two-level systems by Landau |18j . 
Zener [HI] > Stueckelberg [20] , and Rosen and Zener [21] . For an overview see Ref . [22] ■ 

In this work we study the quantum dynamics of the two-qubit gate operations 
described above and use Dykhne's semiclassical result to estimate the probability of 
non-adiabatic transition [23] ■ The applicability of Dykhne's formula is analyzed from 
the standpoint of the theory of semiclassical approximations. These semiclassical 
estimates are found to be in excellent agreement with numerical simulations of the 
exact time evolution. Moreover, in a certain limit our model is integrable, allowing 
us to explicitly calculate and interpret the corrections to Dykhne's formula. 

This paper is organized as follows. Section 2 reviews the derivation of an 
effective two-level model. In section 3, we present our main result — the asymptotic 
estimate of double occupancy, which in section 4 is compared with an integrable 
model and a numerical integration of the Schrodinger equation. In the Appendix, we 
construct the scattering matrix for the integrable model, which has three independent 
parameters. 
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2 Mapping to an Effective Two-level System 



For the purpose of studying double occupancy it is practical to examine the dynamics 
of the quantum gate operation in a subspace spanned by singly and doubly occupied 
states. Following Ref. ^0] with only minor changes of notation, we now detail how 
to reduce the description of a system of two coupled quantum dot spin qubits to an 
effective two-level Hamiltonian. The system is described by a Hamiltonian of the 
form 7i = T + C, where C denotes the Coulomb repulsion between the electrons, and 
T = Yh=12 is the single-particle part with 

The single-particle Hamiltonian hi describes electron dynamics confined to the xy- 
plane in a perpendicular magnetic field B. The effective mass m is a material depen- 
dent parameter. The coupling of the dots (which includes tunneling) is modeled by 
a quartic potential 

VW=V{x, V ) = ^(±(f-<?) % + tf) , (2) 

which separates into two harmonic wells of frequency u (one for each dot) in the 
limit a ^> ao, where a is half the distance between the dots and = Jh/muo is the 
effective Bohr radius of a dot. 

Following Burkard et al. jH] we employ the Hund-Mulliken method of molecular 
orbits to describe the low-lying spectrum of our system. This approach concentrates 
on the lowest orbital states in each dot and is an extension of the Heitler-London 
method jS]. The Hund-Mulliken approach accounts for the fact that both electrons 
can, in the presence of a finite tunneling amplitude, explore the entire system of the 
two dots, and therefore adequately includes the possibility of doubly occupied states. 
In the usual symmetric gauge A = B(— y,x,0)/2 the Fock-Darwin ground state of a 
single dot with harmonic confinement centered around r = (±a, 0, 0) reads 




where Ib = yhc/eB is the magnetic length, and the frequency u is given by uj 2 = 

c<Jq + (ujl/2) 2 where ljl = eB/mc is the usual Larmor frequency. From these non- 
orthogonal single-particle states we construct the orthonormalized states \A) and \B) 
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with wave functions 



<^> = vl _ 2A . g+g > +--w->- <«) 

with S being the overlap between the states (JHJ) and g = (1 — y/1 — S 2 )/S. For 
appropriate values of system parameters such as the interdot distance and the external 
magnetic field, the overlap S becomes exponentially small [Bj. In this limit an electron 
in one of the states \A), \B) is predominantly localized around r— (±a, 0, 0). In the 
following we consider this case and use these states as basis states to define qubits, 
i.e. qubits are realized by the spin state of an electron in either orbital \A), or orbital 
\B). 

An appropriate basis set for the six-dimensional two-particle Hilbert space is given 
(using standard notation) by the three spin singlets 

\Si) = ^= (4 T 4i - 4i4 T ) |o> , (6) 

\S 2 ) = ^(cVi + c B T 4|)i0>, (7) 



\s 3 ) 



^ (4 T 4i - 4 T 4i) |o> , (8) 



and the triplet multiplet, 

l^ 1 ) = 4|4Jo>, (9) 
\t°) = 4(4 T 4i + 4i4 T )|o), (io) 



V2 

IT 1 ) = 4 T 4 T |o). (ii) 

As the Hamiltonian conserves spin, the three triplet states are degenerate eigenstates 
(typically we can ignore possible Zeeman splittings jS]) and have the eigenvalue, 

e trip = 2 £l + V. , (12) 

where we have defined 

e x = (A\h x \A) = {B\hx\B) (13) 
and the expectation value of Coulomb energy 

V- = (T a \C\T a ) , V+ = (Si\C\Si) . (14) 
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An important further observation is that, as a consequence of inversion symmetry 
along the axis connecting the dots, the Hamiltonian does not have any non-zero 
matrix elements between the singlet state | -S3) and other states. Hence, | -S3) is, 
independently of the system parameters, an eigenstate. The eigenvalues of the triplet 
and I S3) , however, do depend on system parameters. The Hamiltonian acting on the 
remaining space spanned by | Si) and | S2) can be written as 

H^ 1 + \ UlI + V + - U -f (I I") (IB) 

where 

t H = -^-[{A\h l \B)+ l -{S 2 \C\S 1 )) (16) 

and 

U H = {S 2 \C\S 2 )-V + . (17) 

The nontrivial part of (J15)) is a simple Hubbard Hamiltonian on two sites and can 
be identified as the Hamiltonian of a pseudospin- 1 object in a pseudomagnetic field 
having a component Uh in the i-direction and Un^H in the x-direction of pseudospin 
space. [Note that this pseudospin is not related to the spin degree of freedom which 
constitutes the qubit!] The basis states themselves are eigenstates only in the case of 
vanishing tunneling amplitude t H where | Si) is the ground state and \S 2 ) is a higher 
lying state due of the Coulomb (Hubbard) energy. In all other cases, the ground 
state has an admixture of doubly occupied states contained in \S 2 ). The energy gap 
between the triplet and the singlet ground state is 

strip - e gs = V„-V + -^ + ^ V / l+^ • ( 18 ) 

A key challenge for state-of-the-art quantum information processing is the con- 
struction of systems composed of two coupled quantum dots which can be coupled to 
perform swap operations lAswi i- e - unitary two-qubit operations which interchange 

the spin states (qubits) of the electrons on the two dots. By combining the "square 
1/2 

root" U sw of such a swap with other isolated-qubit manipulations one can construct 
a quantum XOR gate. A quantum XOR gate, along with isolated-qubit operations, 
has been shown to be sufficient for the implementation of any quantum algorithm 
|25j . Hence a practical and reliable realization of a swap gate would be an important 
step toward the fabrication of a solid state quantum computer. A swap operation in 
the present system is a unitary transformation which turns a state having the qubits 
in different states, say, 

4 T 4J0) = i (|T°> + \S X )) , (19) 
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into a state where the contents of the qubits are interchanged, 



c^\0) = ^= (\T°) - IS,)) . (20) 



V2 



These two states are eigenstates in the case V+ = V- and tu = for which the 
singlet-triplet splitting vanishes. 

As discussed in references El HO], swapping may be achieved by the action 
of a gate that lowers the potential barrier between the quantum dots. This leads to 
exponentially larger values for both V + — V- and tu- It is adequate for our purposes to 
consider a model where V + = V- (consistent with the above limit of small overlap S), 
and the singlet-triplet splitting results entirely from tu- If the duration and amplitude 
of a tunneling pulse is adjusted appropriately, the relative dynamical phase between 
the singlet and the triplet state accumulates a shift of n, 

/oo 
dt (etri P (t) - e ff (t)) = 7T (21) 
-oo 

and the swapping operation between states (pHJj) and (j20j) is performed. However, 
during the operation the state \Si) is coupled to IS2), and they evolve according to 
fTKJ) . Double occupancy errors are thus generically introduced. 

The reduction of the dynamics to the time evolution of a two-level system relies 
on the fact that the system has inversion symmetry along the £-axis in real space 
connecting the dots. This symmetry can be broken if odd powers of the particle 
coordinates Xi are added to the Hamiltonian (£[}, for example the potential of a ho- 
mogeneous electric field. The breaking of inversion symmetry introduces additional 
matrix elements between IS3) and the other two singlets leading to an effective three 
level Hamiltonian. However, as it was shown in Ref. [TO] , this more inclusive Hamil- 
tonian has qualitatively the same properties concerning non-adiabatic dynamics as 
the two-level system on which we shall concentrate in the following. 

So far we have not considered a possible Zeeman coupling to the electron spin. This 
would not change the situation essentially since all states involved in the swapping 
process (|T°), \S\), IS2), and possibly IS3)) have the total spin quantum number 
S z = 0. 



3 Analysis of Non-adiabatic Transitions 

In this section we use Dykhne's formula for non-adiabatic transitions to derive an 
asymptotic expression for the probability of final double occupancy, given physically 
motivated properties of the two-qubit operation. 



7 



As described in the previous section, the modulation of the tunneling barrier 
during the swapping process induces a coupling between the singly occupied qubit 
state l^i) and the doubly occupied state IS2). Their dynamics are governed by the 
effective Hamiltonian 




(22) 

in the 1 5*1,2) basis. The terms omitted from (JT5|l do not contribute to transitions, 
because the identity operator in the \Si^) basis commutes with the remainder of the 
Hamiltonian. The large energy offset Uh between singly and doubly occupied states, 
primarily due to the Coulomb repulsion, is perturbed only by an exponentially small 
additive quantity (proportional to the overlap, S) during the swapping operation and 
is hereafter assumed to be a constant. Our specification of the pulse (Figure ^) 

1 ' cosh(T/2r) 

with dimensionless strength 5 is considered to realistically reflect the tunneling ampli- 
tude that would arise from a modulation of the gate potential [10] ■ The exponential 
dependence of the ramping near t = ±T/2 has its origin in the exponential sensi- 
tivity of the coupling to the gate voltage and in turn the exponential decay of the 
single-particle wavefunctions © in the inter-dot region [20] . The pulse mimics a step 
of duration T and magnitude 5Uh/2, whose ramping on and off have a characteristic 
time r. The perturbation of the instantaneous eigenvalues by the pulse is shown in 
Figure El 

The Schrodinger equation is 

i^ \m) = n cS (t) (24) 

Our task is to find the component of double occupancy in the final state, (^^(oo)), 
given that the prepared state is purely singly occupied, \(Si\ip(—oo))\ = 1. 

Our model involves three dimensionless scales, assigned for our purposes as follows: 
5, X = Uht/2%, and 77 — T jr. Presently, the case of interest is 

A > 1 and 77 > 1. (25) 

The first of these conditions reflects the adiabaticity of the problem. The second 
requires that the ramping on and ramping off of the pulse be temporally well-separated 
and distinct events. 
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Figure 1: A realistic profile of the tunneling pulse, labeled with the characteristic 
time scales. Corresponding to A = 2 and 8 = |, the time scales shown are r and 
T ss 13r. 
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Figure 2: A profile of the eigenvalues ±e(t) corresponding to the pulse of Figure [TJ 



Let us pause and for this paragraph review the familiar notions of transitions 
under the action of a time dependent perturbation. The pulse acts as a transient 
perturbation and otherwise the Hamiltonian (J22|) is diagonal. By force of the adiabatic 
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theorem, the probability of transition among eigenstates vanishes in the limit r — > oo, 
where the ramping on and off of the pulse is adiabatic. In the zeroth order of adiabatic 
perturbation theory, there are no transitions, and the leading behavior of the general 
solution is simply the dynamical phase of each component eigenstate 

mt)) « e^-*^ |a(*)><£i(-oo)|V(-oo)) 

+ e" *£- dt ' e(t,) |6(t)}(6(-oo)|V(-oo)}, (26) 

where |Ci,2 (*) ) are the instantaneous eigenstates (given explicitly in (JMJ)) of Hamilto- 
nian (|2*2*|) corresponding to eigenvalues 

Te(t) = T^Jl + n; , (27) 

respectively. In general, the zeroth order approximation should also include a factor 
representing Berry's phase. However, for a real symmetric Hamiltonian such as ()22)1 . 
Berry's phase is irrelevant, because the Hamiltonian has an inherent planarity. In 
pseudospin one-half notation, 7i cff — H ■ a, the time evolution of the pseudomagnetic 
field H = H(t) is in a plane. If the azimuthal axis (north pole) is chosen to lie within 
that plane, the solid angle subtended by the pseudomagnetic field vanishes identically. 
Although Berry's phase is out of consideration, there are interesting circumstances 
where Berry's phase is relevant to transitions. It can correct the transition amplitude 
pTTj and produce topological selection rules for spin tunneling [321 EH] • O ur problem 
is one of a class initiated by the work of Landau, Zener, and Stueckelberg (LSZ) ^SJ 
UniEOj- However, we emphasize that for our model (with the pulse specified as (|2~3j) ) 
the linearization of Hamiltonian matrix elements near the times where adiabaticity is 
most severely violated is not applicable and leads to an incorrect result. As we will 
see the shape of the pulse is important. 



3.1 Application of Dykhne's Formula 

Returning to our model, we observe that if the time interval t e (— oo, oo) is divided 
into two domains t < and t > 0, and in the limit r\ = ^ ^> 1, the pulse (J23j) is 
approximated by 



t H {t) 



l+e" 



T 
'27 



1+e'i 



T 
'17 



t < 
t > 0. 



(28) 



In each domain the pulse behaves as a step, and the dynamics are integrable (see 
Section 4.) We will focus first on the interval t < 0, where the probability of transition 



10 



to a doubly occupied state P< may be estimated with Dykhne's formula [23] 



|(S 2 |^(0))| 2 ~ e dz£{z \ 



(29) 



where the approximation (|28|) is used implicitly for the instantaneous eigenenergies 
=Fe(t) defined above by (f2~7j). The turning point t — t\, given explicitly below, is 
a complex root of the function e = e(t); in other words, it is an intersection of the 
energy surfaces (curves) of the two instantaneous ("frozen") eigenstates. Our model 
is the patching together of two domains of time, and transitions that occur during 
t < and t > interfere. The expression for the probability of transition during the 
time evolution from t = — oo to £ = oo is 



l(^(oo)>| ; 



i- f, dz e(z) t i , 



i L dz £ ( z ) 



Lxf Ca dze(z) -|9/J ti) dz £ (z) 



i5R f dz e(z) 

e c b e 



Ll 2 ) dz e{z ^ 



(30) 



4 sin 



1 f*2 

- 3« / dz e(z 



P. 



(31) 



where the contours C a> 6 are shown in Figure and according to the sign of the 
integration variable, sgn(3? z), one or the other of the approximations (|28|l is used. 
The turning points t = ti^ appearing in the limits of integration of (|30j) . are chosen 
as the two roots of e — e(t) that are closest to and above the real time axis (see 
Figure EJ) 



T 



i r 



(V — arctan(5)) 



(32) 



They are nonreal because the Hamiltonian ()22j) is nondegenerate for real times. Equa- 
tion (|3"Tj) follows from (jSHJ), because the symmetry of the pulse implies $t(tx) = 3 (£2) 
and P< = P>. The oscillatory first factor of ()31|) is the interference of the dynamical 
phase of each term of (|30p. The magnitude of P is dominated by the second factor, 
P<, whose exponent is given by the following integral 



/■ln(v / T+5 2 ')+i7r— i arctan((5) 

-4A9f / dz 

./ln(^/l+5 ^ ) 



1 + 



5 



l + e z 



-2ttA (l + Vl + S 2 - 5 



(33) 
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Substituting this result in we have 



From (J31)) . we have our main result, an asymptotic estimate for the probability of 
final double-occupancy 

P ~ 4 sin 2 (~ K jf* e(z)j e -^A(i + VI+F-*) ? (35) 

which is shown as a function of 6 in Figure 01 The probability P is characteristically 
nonperturbative in the adiabatic limit r — > oo with fixed, or equivalently A — > oo. 
Hence, the dimensionless quantity associated with the exponential suppression is A 
and has been called the "adiabaticity parameter." For r] = — 1, the approxima- 
tion (|28|) allows us to estimate the argument of the prefactor of (|35|) to exponential 
accuracy, 

- ^J" dz e(z) = VTTP A r] - 2A { In (VlT^ + l) 

- VTTP In (2(1 + 5 2 )) + 5 In [VTT^ + S) 

+ (VT+P - l) ln(5) } + O (e-"/ 2 ) . (36) 

The oscillation with respect to the duration of the pulse T is reminiscent of a similar 
factor in the Rosen-Zener model. The phenomenon of pulsed perturbations that 
return the full amplitude/occupation to the initial state have been studied in the 
context of atom-laser interactions [23 EHl 1201 EOj- In Figures El IH and El we 
compare our semiclassical estimate (|3~5j) with results from numerical simulations of 
the exact quantum mechanical time evolution, following Ref. ^Uj- Both results are in 
excellent agreement and differ only at very small 5, i.e. for weak pulses. Of course, the 
non- adiabatic transition probability vanishes in this limit, whereas the semiclassical 
approximation breaks down (see Section 3.3). This regime is beyond the exponential 
accuracy of Dykhne's formula. The integrability of our model allows us to make 
precise statements about the form and magnitude of the corrections to Dykhne's 
formula (see Section 4). For example, in the limit A ^> 1 and 5X 1 we have from 
the expansion fl63l) that P< ~ (2tt5A) 2 e~ 47rX , while in the same limit the result of 
Dykhne's formula ()35j) gives only the exponential factor e~ AnX without information 
about the prefactor. This explains a trend among Figures El 13 and El namely the 
increasing range, in terms of 8, of validity of Dykhne's formula with increasing A. The 
value for the adiabaticity parameter A = 2, represented in Figure El corresponds to a 
ramp time r = Ah/Un, which was identified in Ref. JTU] as a practical lower bound 
to ensure sufficient adiabatic behavior in a gate operation between two quantum- 
dot spin-qubits. It is interesting that Dykhne's formula remains accurate for smaller 
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Figure 3: The probability for non-adiabatic transitions for A = 2 and rj = 50 as a 
function of 5. We compare our semiclassical estimate according to expression (J35|) with 
results from numerical simulations of the exact quantum mechanical time evolution 
as done in Ref. ^U]. The results are in excellent agreement. 



values of A, in particular A = 1 as seen in Figure El The reason is that the results (|35|) 
and (JoTj) have an incidental factor of 2% in the exponent, giving in practical terms 
the requirement for exponential suppression 27rA ^> 1. 

The expressions ()34jl and (|35|). along with FiguresEllHElcomprise our main results. 
For the remainder of this section we will address the justification and limitations of 
these results. 

3.2 Origin of Dykhne's Formula 

Dykhne derived a concise expression for non-adiabatic transitions from a local analysis 
of the Schrodinger equation in the vicinity of the turning point Dykhne's formula 
can be viewed as a semiclassical approximation, and an elegant interpretation and 
proof was given by Hwang and Pechukas [21] (see also [M].) We will briefly discuss 
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Figure 4: The probability for non-adiabatic transitions for A = 4 and r] = 50 as a 
function of <5. 



the key elements and scope of the proof. Their method was to study the solution 
of the Schrodinger equation in the complex plane of the independent variable, time. 
According to the adiabatic theorem, the projection of the solution onto any eigenstate 
other than the initial eigenstate approaches zero in the adiabatic limit. One might 
suppose that weak statement is all the adiabatic theorem can tell us about transitions; 
however, it does not exhaust its capacities. The reason lies in the following: a basis 
of eigenstates |£i,2(£)), when extended into the complex time plane, is multi-valued. 
In particular, as a basis state is analytically continued across a branch cut of the 
function e = e(t), its long time asymptotics are discontinuously changed. In accord 
with our above two-level problem, we uniquely specify the basis by its asymptotics 

16,2(0)^1^1,2) as t^±oo. (37) 

The multi-valued nature is not manifest on the real time axis, because owing to 
the nondegeneracy of the spectrum ±e(t), the branch points are nonreal. We can 
choose a single-valued basis £1,2^)), which makes reference to |£i,2(0) but has fixed 
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Figure 5: The probability for non-adiabatic transitions for A = 1 and rj = 50 as a 
function of 5. 



asymptotics, by defining rules for continuing the basis states across branch cuts. 
Equivalently, this new basis is said to be defined over a Riemann surface with two 
sheets (copies of the complex time plane), one corresponding to each of the two 
branches of the function \e(t) 2 ] 1 ^ 2 . Crossing a branch cut means passing to the other 
sheet of the Riemann surface. We assign the following relations among the eigenstates 




I6, 2 (t)> 



e 1Q1 - 2 |6,i(*)> 



t G sheet 1 and 
t G sheet 2 



(38) 



where a ij2 are phase conventions that are chosen to maintain continuity of the ba- 



sis 



£i )2 y across the branch cut (no physical quantity will depend on ai j2 -) Given 
— oo)) = oo)), the conclusion of the adiabatic theorem may be restated on a 



Riemann surface as 



£i(t)W) 



vt 



as 



oo, 



(39) 
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Figure 6: The analytic structure of the function e — e(t) shown only in a segment of 
the upper half plane. The contour C a is associated with transitions that occur due 
to the ramping on of the pulse, while contour Cb is associated with the ramping off. 
By Cauchy's theorem, an integral on the contour C a is equal to the integral on the 
contour C' a . Bold lines represent branch cuts, dots represent branch points, and poles 
are denoted with a x . 



where r is the characteristic time scale for variation of Tt e s(t)- The only exception 
to is f° r times within r(9(A~ 2 / 3 ) of a turning point, for there the semiclassical 
criterion pp")) is invalid. As remarked above, the zeroth order approximation (flZo^ 
of the solution as r — > oo is the dynamical phase. The zeroth order approximation 
may be extended into the complex plane by evaluating the dynamical phase on a 
contour C. Continuing with the above example, a state that is purely singly occupied 
at t = — oo is for complex time given by 

&(*)), (40) 

where C is a contour from z = — oo to z — t. The amplitude of transition is readily 
obtained as the projection of the solution, as t —>■ oo on the second Riemann sheet 
(see Figure |7|), onto the doubly occupied state IS2), i.e. 

(S 3 \1>(oo)) = (&(oo)hKoo)) w e 1 ^ e-^C-M), (41) 

where the contour C crosses the branch cut emanating from the branch point that 
is closest to the real axis. Dykhne's formula is simply the square modulus of this 
amplitude. 

In the adiabatic regime, in contrast to the perturbative regime, the leading contri- 
bution to transitions comes from the zeroth order term of perturbation theory instead 



\i/j(t)) w e _i Jc 
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Oi 
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Re(t) 



Figure 7: An example of a Riemann surface with two sheets, a branch point at t — to, 
and a contour C corresponding to a transition. 



of the first order term. By retaining only the zeroth order term, it appears that we 
have neglected completely the coupling t H between states. However, the coupling 
enters implicitly in the multivalued function [e(t) 2 ] 1 ^ 2 and influences the location of 
the turning points — the complex roots of s(t). Transition amplitudes are obtained 
by carefully considering the different branches of this function. In the next section, 
we consider the validity of keeping only the zeroth order term. 

3.3 Validity and Accuracy of Dykhne's Formula 

The theory of semiclassical approximations, especially WKB analysis, provides a foun- 
dation from which to evaluate the validity of Dykhne's formula. The calculation of 
non-adiabatic transitions is closely related to the semiclassical approximation |24j . 
because the semiclassical limit h — > can be mathematically equivalent to the adia- 
batic limit r — > oo. An essential element of the proof by Hwang and Pechukas is the 
existence of a complex time contour that 1) connects the two sheets of the Riemann 
surface and 2) on which the zeroth order approximation of adiabatic perturbation 
theory is the correct leading behavior of the solution in the adiabatic limit. These 
are sufficient conditions for Dykhne's formula to give the correct asymptotic form 
of the transition probability in the adiabatic limit A — > oo. Having established the 
existence of such a contour, one can calculate a more precise value for the prefactor 
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of Dykhne's formula by applying time dependent perturbation theory along the con- 
tour. We expect Dykhne's formula to breakdown when the contour ceases to exist. 
At the limit of its range of validity, the higher order terms become comparable to the 
zeroth order term. Introducing the unitary transformation U that diagonalizes the 
Hamiltonian, i.e. U^HU = £03, we can write the Schrodinger equation in the basis 
of instantaneous eigenstates 

ih± \m = (<t)a s +M(t))m)) (42) 

with the off-diagonal perturbation a(t) = —U* idt U. A dominancy balance among 
the terms gives the condition for the accuracy of the zeroth order approximation 

> \a(t)\ ( 43 ) 



or in scaled time x = t/r, 



h 



»|cs(a;)|. (44) 



h 



For our model of the dynamics, A oc r is the largest scale and |d| ~ 1. The condition 
(JUJ) must be maintained at all points on the contour. Applying (jUJ) on the real axis, 
where \e(x)\ ~ Uh, gives the adiabaticity condition A ^> 1. Additionally, in order to 
connect two Riemann sheets, the contour must pass between two turning points (see 
FigureEI), where \e(x)\ ~ SUh, giving the condition 5X ^> 1. 

Beginning instead from an intuitive approach, we can evaluate the adiabaticity of 
the dynamics along a given contour. To test whether a given contour is adequate, we 
can exploit the analogy between quasi-adiabatic dynamics and semiclassical scatter- 
ing. Recall the semiclassical criterion 

A ( A ) dA 1 
!\ ~ ~dt (45) 

The analog of the de Broglie wavelength X(x) = 2nh/p(x) in scattering problems is 
the period A(t) = 2irh/e(t). In other words, the condition (|45|) says that the change 
of period over the course of one period is small. We now require that the semiclassical 
criterion be obeyed everywhere along an admissible contour, i.e. one that connects the 
two Riemann sheets. To find an admissible contour, we must appeal to the analytic 
structure of the eigenenergy e = s(t), see Figure El For clarity we will focus on the 
time interval t < and operate under the approximation (|28|). The singularities of 
e(t) are branch points at t — — ? — r ln(— 1 ± iS) and poles at t — — ? — r ln(— 1). 
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If we agree to define a branch cut connecting the nearest and next nearest branch 
points to the real time axis 



T 

t\ = — — — t In \A + 5 2 + ir (tt — arctan(<5)) and 

T 



t 3 = — — — r In vl + 5 2 + ir (-7T + arctan(5)) , 

respectively, then an admissible contour is one that crosses this branch cut exactly 
once. For the semiclassical criterion to be obeyed, the admissible contour cannot 
pass too close to a branch point. In essence, if 6 is too small, the contour is pinched 
between the branch points t\ and t 3 . Evaluating the maximum of dA/dt on a contour 
C that crosses the branch cut between t — t\ and t = t 3 we arrive at the condition 
5 ^> A -1 . Together with the adiabatic limit A ^> 1, we have the following conditions 
on the interdependence of the physical parameters 

Uttt 

A J- > 1 and (46) 

n 

SX^^l > i. (47) 

Each of these dimensionless quantities is a product of a characteristic energy and 
time scale. If these conditions are not satisfied, there does not exist a contour on 
which the motion is adiabatic. The integrability (Section 4) of our model allows us 
to investigate the intermediate regime A ^> 1 and SX 1, where Dykhne's formula 
cannot be justified with the analysis of Hwang and Pechukas. 



4 Identification with an Integrable Model 

The result obtained by Dykhne's formula in Section 3.1 is now shown to be equivalent 
to the exact result for an integrable model in the appropriate limit. 

Under the approximations (J28)) for the time intervals t < and t > 0, the Hamil- 
tonian 



is approximated by 



H eS « I < 49 
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with 



n< = -^f 1+e lV" ] (so) 



2 



1+e t 27 



The Hamiltonians W< and 7^> can be obtained as a special case of 

-j ( b a + c tanh|\ . , 

^ exact ~ ^ a + c tanh § -6 J ^ 

by identifying ±c = a = — 5A/2, 6 = —A and rescaling time x = t/r ± T/2r, respec- 
tively. The Schrodinger equation 

ld x \${X)) = Hexact |V0&)> (53) 

is exactly solvable jSH] (see Appendix.) 

In analogy with one- dimensional scattering, the transition amplitude from a singly 
occupied state \Si) to a doubly occupied state \S 2 ) may be viewed as an off-diagonal 
element of the scattering matrix S that connects the coefficients of the asymptotic 
final states to the asymptotic initial states. The asymptotic states are the limit as 
t — > ±00 of the instantaneous eigenstates |£i,2(*)) of H e fi corresponding to eigenvalues 
=Fe(t), respectively, 



™> -M*3Y (54) 



The leading behavior of the long time asymptotics of a general solution has the form 
( ai e -</' dt ' ( - £(t,)) |^i (t)> + a 2 e dt ' ^ |f 2 (t)) as t -> -00 

\m) = (55) 

1 h e - l f dt ' |a(t)> + &2 e dt ' £(t,) as t -> 00 

and the scattering matrix relates the coefficients 

(£)-*(:)■ « 56) 
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Referring to the statement of our problem in Section 3, the amplitude of final double 
occupancy is the element S21 of the scattering matrix, which we parametrize as 



S 



— °12 °11 



(57) 



The scattering matrices S < and S > associated with the Hamiltonians 7i < and 7i > 
may be obtained by substitution from the exact scattering matrix W (derived in the 
Appendix) associated with the Hamiltonian 7i ex act- By the symmetry of the pulse, 
we have S > = 5<. Patching together the two domains of time evolution (pf9~|) . we find 
the scattering matrix 



S = S{ exp ^ i 0-3 $i J* 2 dt e(t)J <S<, 



(58) 



where the integral of the exponent has been estimated in (j3T?j) . and the elements of 
S < are obtained from W 



P<)n 



r(i2A) r(i2/i) 



(s. 



O12 



I 2/i 

[i + A r(i/i + iA + i<5A) T(i/i + iA - \b\) ' 

/ 2/i r(i2A) r(-i2/i) 

/i - A r(— iju + iA + \b\) T(-i/i + iA - i<5A) 



(59) 



(60) 



where fi = AVi + 5 2 . Dykhne's formula for P K is recovered as exactly the leading 

|2 



term of | (S. 



<J2l\ 



in the limit A, 5X ^> 1, 

|2 



1(5. 



<J21\ 



sinh 


ttA 


yi + 5 2 - 1 + 5) 


sinh 


-ttA 


[V 1 + 5 2 - 1 - 5) 


sinh (27rA) sinh 


2ttAV1 + S 2 j 





27rA(l+v / I+F-<5) _ e -27rA(<5+l-VI+F) + \ 



(61) 



(62) 



The nonperturbative corrections are typically very small. For our problem A = 2, 
8 = 1/2, and the relative contribution of the second term of ()62|) is less than one 
percent. This accounts for the excellent agreement in Figures and0 between 
the probability of double occupancy as given by the semiclassical result (j3*B]) based 
on Dykhne's formula and the result of a numerical integration of the Schrodinger 
equation. We can interpret the subleading term of (J62|) as the contribution from the 
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Figure 8: The contour C gives a subdominant correction to Dykhne's formula. 



contour C of Figure |H1 which crosses the branch cut three times. The sign of the 
correction is negative and arises from the factor e iai e ia2 associated with matching 
the basis (|3*8*|) across the branch cut. For our problem ot\ = and = 7T. Similarly, 
it may be possible to extend Dykhne's formula and reproduce the exact result (|HT]) by 
summing over all homology classes of the Riemann surface and within each homology 
class over all paths that give distinct nonzero values for 9 J c dt e(t) [35]. For many 
physical problems this type of nonperturbative correction is dominated by perturba- 
tive corrections along the contour — those mentioned in Section 3.3. The striking 
absence of perturbative corrections in the limit A, 5X ^> 1 is a unique artifact of the 
integrability of our model. 

With a knowledge of the exact result, we can also investigate the intermediate 
regime A > 1 and 5X <C 1, where the analysis of Hwang and Pechukas (Section 3.2) 
cannot be used to prove Dykhne's formula. In this limit, the transition probability 
(pjTj) becomes 

P< ~ (2vr5A) 2 e~ 47rA . (63) 

Although Dykhne's formula does not apply in this limit, because <5 _1 and not A is 
the largest scale, it nevertheless gives the correct controlling factor e~ 47rA of (|63|). 
except for 5 that are exponentially small with respect to A. This exponential factor 
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is resilient and remains the controlling factor for a range of parameters beyond the 
naive expectation based on the arguments of Section 3.3. 

5 Conclusions 

The dynamics of two coupled quantum dot spin-qubits can be mapped to an effective 
two-level system, where non-adiabatic transitions correspond to double occupancy. 
We have estimated the probability of final double occupancy with Dykhne's formula. 
In the adiabatic regime, the pervasive feature of transitions is their exponential sup- 
pression by a dimensionless adiabaticity parameter, A. Our main result (J33j) was 
expressed in terms of the dimensionless quantities A, 5 and r\. An integral constraint 
(12 1|) on the swapping operation gives one relation among the three dimensionless pa- 
rameters. The problem is uniquely defined by specifying any two, and in a solid state 
setting, conservative estimates are A ~ 2 and 5 ~ 1/2. The probability of double 
occupancy P « 10 -10 is sufficiently rare that the operation of a quantum gate will 
not be obstructed by this type of error. It is noteworthy that the probability of double 
occupancy (J3T)|) has nodes for 

krr = ft [' dt e(t) « Vl + 5 2 Xr l k e Z. (64) 

Jti 

However, this property is not immediately relevant to the suppression of transitions, 
because the oscillatory factor, sin 2 (Vl + ^Xr/^j , of (|3*5|) vanishes algebraically and 
for it to provide an improvement upon the exponentially small factor from Dykhne's 
formula, the argument would have to be tuned exponentially close to kn. Thus, 
naively the errors associated with inaccuracies in satisfying the integral swapping 
constraint (}2~Tj) will be much greater than double occupancy errors. Other important 
sources of error are dephasing and decoherence of the qubit states. 

We have reviewed a physically motivated derivation of Dykhne's formula [21]. The 
theory of semiclassical approximations underlies Dykhne's formula and its validity is 
appropriately judged within that framework. The semiclassical estimates obtained 
from this approach are in excellent agreement with numerical simulations of the full 
quantum mechanical time evolution. The corrections to Dykhne's formula are of two 
types: perturbative and nonperturbative. The former appear to vanish for integrable 
models, and we have interpreted a nonperturbative correction as the contribution of 
a contour in the complex time plane [3*5] . 
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A Integrable Dynamics 

The following Hamiltonian which possesses three independent parameters 

<U ( \ l b a + C tann f \ 

1^ + ctanhi _ 6 2 J (65) 

has integrable dynamics in the sense that the Schrodinger equation, 

(i d x - ftexactOr)) |V(a;)) = with \^{x)) = ( ^ ) (66) 

has a solution in terms of special functions. Our aim is to construct the scattering 
matrix W associated with the dynamical problem. For clarity of presentation, we will 
consider the Hamiltonian 

/ a + c tanh § b \ 

"W = ( 6 2 -(a + ctanh f ))> < 67 > 

which differs from 7i ex act by a constant unitary transformation V 

V ] Kxact V = W With V =±=(\ \ y (68) 

Introducing the mixing angle <fr defined as 

^ a + c tanh | ^ 
the instantaneous eigenstates |xi,20e)) corresponding to instantaneous eigenvalues 



X N 2 



±e(x) = ±\jb 2 + (a + ctanh-) (70) 
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respectively, are parametrized as 



cos • 



1 / Je + (a + ctanh f) . . vw , . 

Xi (*)) = -7= V / - — = I and 71 

IA W/ v 7 ^ I - (a + c tanh f) J V sm f 



i*w> - i y ^+^T ] - ( ~ sn > i • (72) 



2e \ Je + (a + ctanh |) y \ cos I 

The two-state Schrodinger equation may be converted to a second order differential 
equation for the components Ci >2 = Ci j2 (o;), 

^ci )2 + Qi, 2 (x) ci, a = (73) 

where 



Qi,2(x) = 6 2 + [a + ctanh - J ± i - c seclT-. (74) 

Changing the dependent and independent variables as c\p(x) = z ±lv (z — l) 1M wi )2 (-2) 
and z = \ (l + tanh |j, transforms (J75J) into the standard form of the Gauss hyper- 
geometric equation, for example, for wi(x) 

z(l — z) d 2 z w\ + [k — z(i + j — 1)] d z wi — ij w% = 0, (75) 

where the arguments and exponents are defined as follows 

% = i// + \v — i2c 
j = ifx + iu + i2c + 1 
k = i2z/ + l 



v = e(—oo) = \/b 2 + (a — c) 2 



fj, = e(oo) = yb 2 + (a + c) . 

The two linearly independent solutions Ui(z), Vi(z) of the differential equation for 

101(2:) are 

Ut(i,j,k;z) = z w (z-l) lfl 2 Fi(i,j, k;z) and (76) 

vi(l, m, n; z) = z x ~ k U\{% — k + 1, j — k + 1, 2 — k; z) 

= z~' w (z-lf* 2 F 1 (i-k + l,j-k + l,2-k;z), (77) 
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where 2 Fi k; z) is the Gauss hypergeometric function [37 . 

The amplitude of transition may be viewed as the off-diagonal element of the 
scattering matrix W that connects the asymptotic final and initial states for t — ► ±00 
(see ()56J1 .) The scattering matrix is parametrized as 



W u W12 



(78) 



because it has the properties W = 1 and det(W) = 1. From the asymptotics of 
Ui, «2> Vi, and t>2 in the limits z — > and z — > 1, we find 



11 



r(i - \2v) r(-i2//) 



& 2 + (^-(a + c)r 

\J 6 2 + (1/ - (a - c)) 2 r (1 - i (// + v - 2c)) r (-i fa + v + 2c)) 



(79) 



12 



& 2 + (/i+(a + c)y 



T(l - \2v) T(i2fi) 



\b 2 + 



v — (a 



c)y V(i{n-u- 2c)) r (1 + i (fi - v + 2c)) 



(80) 



and 



|W- 



11 



12 



sinh 7r (2c + + v) sinh 7r (—2c + fi + z/) 
sinh 27r/i sinh 2ttv 

sinh 7r (2c — fi + v) sinh 7r (2c + fi — z/) 
sinh 27r/i sinh 2nu 



(81) 
(82) 
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